Aging Effects Across the Metal-Insulator Transition in Two Dimensions 
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Aging effects in the relaxations of conductivity of a two-dimensional electron system in Si have 
been studied as a function of carrier density. They reveal an abrupt change in the nature of the 
glassy phase at the metal-insulator transition (MIT): (a) while full aging is observed in the insulating 
regime, there are significant departures from full aging on the metallic side of the MIT, before the 
glassy phase disappears completely at a higher density n g ; (b) the amplitude of the relaxations 
peaks just below the MIT, and it is strongly suppressed in the insulating phase. Other aspects of 
aging, including large non-Gaussian noise and similarities to spin glasses, also have been discussed. 

PACS numbers: 71.55.Jv, 71.30.+h, 71.27. +a 



The glassy freezing of electrons, resulting from 
the competition between disorder and strong electron- 
electron interactions, may be crucial for our understand- 
ing of the behavior of many materials near the metal- 
insulator transition (MIT) [1]. The dynamics of such 
electron or Coulomb glasses [2] , however .remains poorly 
understood, and experimental studies 0,0, H[ are still rel- 
atively scarce. In two dimensions (2D), moreover, even 
the very existence of the metal and the MIT have been 
questioned. Although there is now considerable experi- 
mental evidence in favor of such a transition [6] , there is 
still no general agreement even on the fundamental driv- 
ing force for the 2D MIT. Theoretical proposals describ- 
ing the 2D MIT as the melting of a Coulomb glass [7|, |8|, |9[ 
have found support in observations and studies of glassi- 
ness in a 2D electron system (2DES) in Si. 

Aging [14, 15], one of the key characteristics of glassy ma- 
terials, also has been demonstrated in this system [13 ]. 
but its properties have not been investigated in detail. 
Surprisingly, although consistent with theory H, glassi- 
ness sets in on the metallic side of the MIT [To|, Ujl E3] , 
i.e. at an electron density n g > n c (n c - the critical den- 
sity for the MIT), thus giving rise to an intermediate, 
poorly metallic and glassy phase between the metal and 
the (glassy) insulator. Here we focus on aging effects, 
which have been instrumental as a probe of complex 
nonequilibrium dynamics in many types of materials. We 
find an abrupt change in the aging properties that occurs 
at n c . The observed change in the dynamics at the 2D 
MIT itself, in addition to that at n g [T3, [n|, [T3|, might 
help to determine the validity of different theoretical pic- 
tures. 

The system is said to exhibit aging if its response to 
an external excitation depends on the system history in 
addition to the time t. In electron glasses, aging has been 
studied mostly and most easily [4( by looking at the relax- 
ations of conductivity a(t) towards its equilibrium value 
do after a temporary change of electron density n s dur- 
ing the waiting time t w . Aging is observed if t w <C r eq 
(r eq - equilibration time), i.e. if the system is not able 
to reach equilibrium under the new conditions during t w . 



It is manifested in the dependence of a(t) on t w such 
that, in those strongly localized systems, the aging func- 
tion a(t,t w ) is just a function of t/t w [4]. This is known 
as simple, or full aging. It is interesting that, in spin 
glasses, full aging has been demonstrated only relatively 
recently [lf|. In general, however, the existence of a char- 
acteristic time scale t w does not necessarily imply simple 
t/t w scaling [ljj. In the mean- field models, in fact, two 
different cases are distinguished: one, where full aging 
is expected, and the other, where no t/t w scaling is ex- 
pected [l8(| . Experimentally, departures from full aging 
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temperature T — > (hence glass transitio n T„ = 0) [l2[, 
aging was also observed for t w <C r eq (T) [13[ using the 
experimental protocol described above. The goal of this 
study is to investigate that aging regime in detail. In 
particular, unlike previous work [13[, here T is kept fixed 
at 1 K such that r eq is astronomical [19[ and the 2DES is 
always deep in the t w <C r eq limit; <j(t,t w ) are then ex- 
plored systematically both as a function of final n s and of 
the difference in densities during and after t w . Our main 
results include: a) a(t) obey a power-law dependence for 
t <t w and a slower relaxation law for t>t w , thus show- 
ing a memory of the time ("age") t w ; b) a(t,t w ) exhibit 
full aging for n s < n c ; c) as n s increases above n c , there 
is an increasingly strong departure from full aging that 
reaches maximum at n s — n g ; d) for a given tqjj , the am- 
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interesting and surprising suppression of the relaxations 
on the insulating side of the 2D MIT. 

The experiment was performed in a 3 He system 
(base T = 0.24 K) on the same (lOO)-Si metal-oxide- 
semiconductor field-effect transistors that were used in 
previous studies [HI, [l3[. The two devices (A and B, 
with 1 x 90 /im 2 and 2 x 50 fim 2 length x width, re- 
spectively) had a 50 nm oxide thickness, and a rela- 
tively large amount of disorder (the 4.2 K peak mobility 
w 0.06 m 2 /Vs with the substrate (back-gate) bias [20[ 
V su b = —2 V). n s was varied by the gate voltage V g , such 
that n s (10 n cm- 2 ) = 4.31(V g [V] - 6.3); n^(10 n cm- 2 ) = 
(7.5 ± 0.3) and n c (10 n cm- 2 ) = (4.5 ± 0.4), where n g 
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FIG. 1: (color online) (a) V g (t) and T(t) in a typical ex- 
perimental protocol, which always starts with the 2DES in 
equilibrium at 10 K [13|]- The results do not depend on the 
cooling time (varied from 30 minutes to 10 hours), nor on t± 
and £2 (varied from 5 minutes to 8 hours each), (b) The 
corresponding cr(t). The relaxation of a during t w is too 
small to be seen on this scale, (c) a(t > 0) for several t w , 
as shown; V = 7.0 V [n (10 n cm" 2 ) = 3.02 < n c ], Vi = 11 V 
[ni(10 n cm -2 ) = 20.26]. The data are not shown for the first 
few seconds (comparable to our sampling time ~ 1 s). The 
dotted lines are linear fits for t < t w . (d) The same data as 
in (c) but plotted vs. t/t w . The dotted line is a fit for t < t w 
with the slope —a = —0.081 =b 0.005. Inset: cr(t = 1s)/<jo vs. 
t w . The dotted line is a fit with the slope a = 0.076 =b 0.005. 



was determined from the onset of slow, correlated dy- 
namics in noise and n c from <r(n s ,T) measurements on 
both metallic and insulating sides [IcJ EH • The devices 
and the standard ac lock-in technique (typically 13 Hz; 
5-10 /iV excitation voltage) were described in more de- 
tail elsewhere [10]. The two samples exhibited an almost 
identical behavior. Unless noted otherwise, the data pre- 
sented below were obtained on sample A. 

We employ the so-called "gate protocol" [4], where V g 
is changed from an initial value Vq (density no), where the 
2DES is in equilibrium, to another one, V\ (density 774), 
where the system attempts to equilibrate during t w , but 
tw <^ T eq [Figs. [2(a) and [2(b)] . After V g is changed back 
to Vq at t = 0, the relaxation of the "excess" conductivity 
a(Vo,t)/<j(Vo) is studied. Figure [2(c) shows some a(t) 
measured for several t w and no < n c . Although n\ > n„ : 
small amplitude relaxation during t w is still visible [I2L 
13|, [2l|. It is clear that t w has a significant effect on cr(t). 
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FIG. 2: (color online) Noise in a for several no, before and 
after a change of n s to m during t w . The slowly relaxing 
background < a(t) > has been subtracted from the upper 3 
"after" curves. The data are shifted for clarity. For the two 
lowest no, the signals measured after t w were divided by 10. 



In fact, all the cr(t, t w ) data can be collapsed onto a single 
curve simply by rescaling the time axis by t w [Fig. [2(d)]. 
Therefore, in this case, the system exhibits full aging at 
least up to t « (2-S)t w . We note that the relaxations can 
be described by a power law a(i)/ao ex (t/t w )~ a for times 
up to about t w , followed by a slower relaxation at longer 
t. This means that the memory of t w is imprinted on the 
form of each cr(t). As a consistency check, the fits to the 
individual <j(t) for t < t w [see Fig. [2(c)] yield the initial 
amplitudes of the relaxation a(t = 1 s)/<to cx t% [22] with 
the same value of a (Fig. [2(d) inset). 

The quality of the data collapse is difficult to estimate 
with high precision since a(t) are quite noisy. We ar- 
gue that the large non-Gaussian noise seen in Fig. [2(c) is 
intrinsic to aging, similar to other materials out of equi- 
librium [13, Indeed, Fig. [2] shows that, for n s < n g , 
the noise after a temporary n s change is much larger than 
the noise before the change. At the lowest no, this differ- 
ence amounts to more than an order of magnitude. While 
a more detailed comparison of the "before" and "after" 
noise will be presented elsewhere, we note the following. 
First, the noise in the aging regime grows with decreas- 
ing T, while at the same time the relaxations a(t)/ao 
become smaller. Therefore, in order to optimize the sig- 
nal to intrinsic sample noise for the study of aging, most 
of the data were taken at 1 K. Second, they show that 
the previous noise studies in a 2DES MM 

were actu- 
ally performed in the aging regime, since V g was varied at 
T - 1 K albeit in small steps (AV g =0.1 and 0.01 V, or 
smaller, in Refs. [l(j and [HI ], respectively). Such small 
AV^ did not produce any visible relaxations within the 
noise, consistent with the results below (Fig. SJ). 

Full t/t w scaling is exhibited for all n s < n c . However, 
as soon as n s > n c , we find systematic deviations from 
full aging [e.g. Figs.[3Ka) and [31(c)] . In other glassy mate- 
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FIG. 3: (color online) (a), (c), (e) Relaxations for different 
no(10 n cm -2 ) and a fixed m (10 n cm _2 ), scaled with the wait- 
ing time t w . (b) , (d) Scaling with improves the collapse of 
the data, (f) \i vs. no. \i does not depend on n\. 



rials 14, 24|, it was found that the data could be scaled 
with a modified waiting time {t w )^ , where \i is a fitting 
parameter (/i = 1 for full aging). Even though \i may not 
have a clear physical meaning, the /i-scaling approach has 
proved to be useful for studying departures from full ag- 
ing 0, E3- By adopting a similar method, we find that 
it is possible to achieve an approximate collapse of the 
data [Figs.EKb) and^d)]. The plot of /i vs. n [Fig. 13(f)] 
shows a clear distinction between the full aging regime 
for n s < n c , and the aging regime where significant de- 
partures from full scaling are seen. It is striking that the 
largest departure occurs at n s ~ n g . For n s > n gi it ap- 
pears as if the full aging is restored, but this may be an 
artifact of trying to collapse very small (cr/cro < 2 per- 
cent) relaxations accompanied by instrumental (white) 
noise of comparable magnitude [Fig. [31(e)] . There are no 
relaxations, within the noise, for Vq > 9 V at T = 1 K. 
We have also determined that fi does not depend on T. 

Another important issue that needs to be addressed is 
the role of V\ and of the step size AV g = (Vi — Vq). While 
\i depends only on Vo (provided t w <C r eq at V\ ) , the am- 
plitudes and the slopes of the relaxations for a given t w 
do depend on AV g . Figure [4](a) shows such an example, 
where <j(t) are presented for fixed Vo, t Wl and different V\. 
The data can be fitted by <j(t)/<Jo = [cr(t = 1 s)/<jo] t~ a 
(t in units of s) at least up to t ~ t w . The effect of 
AV g on the slope a and the amplitude a(t = 1 s)/ao is 
given in Figs.H^b) and[4|(c), respectively, for different Vo 
and t w . The increase of both quantities with AV g is not 
too surprising, since larger AV g lead to a higher num- 
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FIG. 4: (color online) (a) a(t) for fixed no(10 n cm -2 ), t w , 
and several m (10 11 cm -2 ). The dashed lines are fits, (b), (c) 
Slopes a and amplitudes cr(t = 1 s)/ao of the relaxations, 
respectively, vs. An s (10 11 cm -2 ) for two fixed no and several 
t w . The values corresponding to no(10 11 cm -2 ) = 6.46 were 
multiplied by 10. Open symbols: device B. Dotted lines guide 
the eye. (d) Sample B. a(t) for small AV g scaled with t w . 



ber of new electrons in the 2DES during t Wl taking the 
system farther away from its original state at Vq. It is 
intriguing, though, that the results [Figs. 2Kb) and[4](c)] 
seem to suggest that, for very small AVg, both a and 
log[cr(l s)/(Jo] may become negative. This corresponds to 
<j(t) approaching <jq from "below" instead of from above 
as in Fig. [He) [25[. Indeed, we have observed instances 
of such behavior, as shown in Fig.[4Kd) for the case where 
both Vo and V\ happen to be in the insulating regime. 
The nonmonotonic dependence of the system response, 
here characterized by a and a(l s)/cr , on the perturba- 
tion AV gi bears a remarkable resemblance to the results 
of negative T cycle spin glass experiments (see Figs. 1-3 
in Ref. [26[), where the so-called "memory anomaly" at 
low AT was interpreted [26] based on the hierarchical 
distribution of states in the free energy landscape. The 
noise studies in the aging regime [11] have already pro- 
vided support for the hierarchical picture of glassiness in 
the 2DES. Unfortunately, it is not possible to perform 
systematic studies of aging for small, including negative, 
AVg in a 2DES, because then a(l s)/ao is typically com- 
parable to the large noise that is intrinsic to aging. 

The raw data [e.g. Figs. C3(a), Etc), andC3(e)], as well 
as Figs. HKb) andHtc), indicate also large variations in 
<j(l s)/(7o and a as a function of Vo for fixed t w and AV g . 
This was explored in detail by varying Vo, while keeping 
t w and V\ fixed [Fig. 03(a)]. We expect, based on Fig. [H 
to see an increase of both quantities with decreasing Vo , 
since AV g is getting larger. While an increase is observed 
[Fig. (5](b)], it is much stronger than expected. For exam- 
ple, even for the same An s (10 n cm -2 ) = 13.79 and t w = 
1000 s, log[cr(no(10 n cm- 2 ) = 3.88, t =1 s)/cr ] « 0.4 
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FIG. 5: (color online) (a) a(t) for fixed t w , n±, and several 
no. (b) Relaxation amplitudes vs. no for several t w . Dotted 
lines guide the eye. The error bars typical of low and high no 
values are also shown. 



[Fig. Sfc)], i.e. still an order of magnitude larger than 
for no(10 11 cm -2 ) = 6.46. Even more striking, however, 
is the dramatic decrease of cr(l s)/<Jo [Fig. 03(b)] as no is 
reduced further, contrary to simple expectations based on 
Fig. [H The slopes a exhibit the same nonmonotonic be- 
havior [Fig. 03(a)]- Again, the abrupt change in the aging 
properties occurs at n s ~ n c : the relaxation amplitude 
peaks just before the system becomes metallic, and it is 
suppressed upon going deeper into the insulator. This is 
reminiscent of the suppression of the fluctuations of a(t) 
observed earlier in the same system at n c (Ref. [10[, Fig. 4 
inset). While a clear understanding of these two effects 
is lacking, it is plausible that collective charge rearrange- 
ments that are responsible for the slow dynamics will be 
suppressed as the 2DES becomes strongly localized. It 
is interesting to speculate whether this is related in any 
way to the problem of many-body localization [27[. It 
would be also of interest to study aging deeper in the in- 
sulator, in the variable-range hopping regime, but that is 
not possible because of the small relaxation amplitudes 
and the large intrinsic sample noise. The effects of disor- 
der could be explored further by extending the relaxation 
studies to cleaner 2DES, where n g > n c 

In summary, a detailed study of aging in a 2DES in Si 
shows an abrupt change in the nature of the glassy phase 
at the 2D MIT before it vanishes entirely at a higher 
density n g . The results put constraints on the theories of 
glassy freezing and its role in the physics of the 2D MIT. 
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